Lithium Diffusion in the Post-Recombination Universe and Spatial Variation of [Li/H] 
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The observed amount of lithium for low metallicity population II stars (known as the Spite 
plateau) is a factor of ~ 3 — 5 lower than the predictions of the standard cosmology. Since the 
observations are limited to the local Universe (halo stars, globular clusters and satellites of the Milky 
Way) it is possible that certain physical processes may have led to the spatial separation of lithium 
and local reduction of [Li/H]. We study the question of lithium diffusion after the cosmological 
recombination in sub- Jeans dark matter haloes, taking into account that more than 95% of lithium 
remains in the singly-ionized state at all times. Large scattering cross sections on the rest of the 
ionized gas leads to strong coupling of lithium to protons and its initial direction of diffusion coincides 
with that of H + . In the rest frame of the neutral gas this leads to the diffusion of H + and Li + out 
of overdensities with the trend of reducing [Li/H] in the minima of gravitational wells relative to 
the primordial value. We quantify this process and argue that, with certain qualifications, it may 
have played a significant role in creating local lithium deficiency within the primordial dark matter 
haloes, comparable to those observed along the Spite plateau. 



The primordial abundances of light elements, 4 He, 
D, and 7 Li, offer unique window into the very early 
Universe at redshifts of z ~ 10 9 . In recent years, 
this probe has been sharpened: the only free param- 
eter that enters the standard big bang nucleosynthesis 
(SBBN) calculations - baryon-to-photon ratio ^ - has 
been measured to great accuracy via the CMB experi- 
ments. The prediction for the primordial fraction of 7 Li 
is [ 7 Li/H] S BBN = (5.07i^2) x 1CT 10 (see e.g. [I]), which 
is a factor of 3 — 5 smaller than the Spite plateau value of 
7 Li, (1.23^0 ^) xlCT 10 0, an observationally determined 
value of the lithium abundance in the atmospheres of hot 
Population II halo stars. It is unclear whether stellar de- 
pletion of Li could account for such a large deficit, and 
speculations of non-standard physics being behind the 
discrepancy flourished (see, e.g. reviews [3]). Most re- 
cently, this lithium problem has been further complicated 
by the observation of the deterioration of the plateau at 
the lowest metallicities, Z < 1.5 x 1CP 3 , where the dis- 
crepancy with SBBN value becomes even larger [3 . This 
might be pointing towards additional "missing" pieces 
of physics unrelated to the stellar physics, starting from 
evolution of primordial gas leading to the formation of 
PopII stars with lowest metallicities. 

It is important to realize that the observations of 
lithium abundance along the Spite plateau reflect the 
"local" formation environments of the oldest stars in 
our Galaxy, while the SBBN predictions are "global". 
Non-standard cosmology with an O(l) downward fluctu- 
ation of baryon-to-photon ratio in the patch of the Uni- 
verse that includes Milky Way can give [ 7 Li/H]i oca i < 
[ 7 Li/H]sBBN [5]. However, the standard physics processes 
may also lead to the local under- or over-abundance 
of lithium relative to the SBBN prediction. In stan- 
dard cosmology, the spatial fluctuations of 7 Li/H are ini- 
tially small, but consequently amplified by the growth of 



structure in combination with diffusional processes in the 
Early and Late Universe. 

The existing astrophysical literature covers the diffu- 
sion of elements in stellar atmospheres and in the clusters 
of galaxies [SHE]- In contrast, the studies of primordial 
element diffusion in the early Universe are very sparse. 
Ref. [9] addresses the evolution of elemental abundance 
in the linear regime, Sp/p <C 1. Although the linear 
regime by definition does not allow for large effects in the 
abundances, [5] find that qualitative trend is such that 
lithium, owing to its larger mass, tends to accumulate 
more in the minima of gravitational potentials compared 
to hydrogen. Since the star formation should also occur 
inside gravitational wells, the qualitative trend inferred 
from [S] is [ 7 Li/H]i oca i > [ 7 L1/H]sbbn, which does not 
help to solve the lithium problem in any way. There is, 
however, an important assumption made about the neu- 
trality of lithium in [9J , which does not hold in the early 
Universe. In fact, after the H recombination, lithium ex- 
ists predominantly in the singly-ionized state, Li + [10] . 
There are two reasons for that: firstly, the 5.39 eV ioniza- 
tion potential for lithium means that the recombination 
temperature is smaller than that of H by the factor of 
~ 2.5, at which time the density of the free electrons is 
depleted, and the recombination rates for Li are less than 
the Hubble expansion rate. In addition, the non-thermal 
population of photons from residual e — p recombination 
causes photo-ionization of neutral Li fraction and keeps 
its abundance below a few percent level throughout the 
cosmic history all the way to reionization at z ~ 10 |10) . 

In this paper, we show that the fact that lithium re- 
mains in the Li + state has direct consequences for its 
diffusion after hydrogen recombination. In particular, 
we show that owing to the large scattering cross section 
on protons, Li + stays spatially bound to H + , and the di- 
rection of their diffusion is against the gravitational force 
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in the rest frame where the neutral hydrogen. From our 
analysis it follows that [ 7 Li/H] grav m ; n < [ 7 Li/H]sBBN, 
which can have significant implications for the cosmolog- 
ical lithium problem. In the rest of the paper, we expand 
on this observation in some detail. 

Direction of Lithium Diffusion. We consider the 
equations for cosmological fluids of different primordial 
species with number densities n a , where a spans H, 4 He, 
e, p, 7 Li. Note that after most of the hydrogen becomes 
neutral, the recombination rate for residual e and p is 
much smaller than the Hubble expansion rate and they 
can be treated as separate species. The presence of small 
quantities of D and 3 He will not affect the evolution of 
7 Li. Thus, we have the system of equations for the aver- 
age velocities V a of individual species, 
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In these equations, g is the gravitational acceleration, E 
is the electric field strength, P ai p a = m a n a , and q a are 
the partial pressure, mass density, and the electric charge 
for different species respectively. For the electromagnetic 
effects we assume the tight charge coupling approxima- 
tion. The V a — diffusion term is governed by the 
diffusion coefficients r 
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that are in turn determined by the transport cross sec- 
tions o a b averaged over the microscopic velocity distribu- 
tion. We use lower and upper case to distinguish between 
thermal v and diffusional V velocities. [i a b is reduced 
mass and T is the temperature of the matter species. 
Finally, the last term in ([I]) accounts for the possibility 
of additional external forces, such as radiation pressure, 
Lorentz force, etc., with dependence on species index a. 
However, we take f° xt = for the rest of this analysis. 

In the next step, we solve Equations ([!]) in the regime of 
small density perturbations, 6p a /p a 5s 1, and specifically 
consider a sub-Jeans regime for baryons, that are forced 
inside an already formed dark matter halo by its gravita- 
tional acceleration g. For a realistic choice of parameters 
I/r — > oo is a good zeroth order approximation, leading 
to vanishing Vj in hydrostatic equilibrium. Thus, assum- 
ing that initial distribution of elements is uniform, one 
gets a relation between gradients of individual pressure 
contributions and g (see, e.g. [7]), 
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and the quasi-static version of Eq. ([I]) reduces to a set 
of algebraic equations, 
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We assume 25% mass fraction of 4 He so that ri4 He /n H = 
1/12 and m~ (4m p +12m p )/(l+12) = jfm p . The reduc- 
tion of m due to ionized fraction can be safely neglected. 

For the two dominant neutral components, hydrogen 
and helium, the solution is readily found: 
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One can see that V Ho — V H is parallel to g, as expected. 

We now turn to the diffusion of charged particles and 
account for E in the equations. Solving them for electrons 
with the use of ^ and m e -C m atom , one can easily find 
E ~ — VP e /(en e ) — — (m/e)g, where e is the positron 
charge. Carrying this to the equation for protons (or 
H + ), we get the solution for the relative diffusion velocity: 
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using Equation (|5j). The appearance of m p /2 in this 
equation is easy to interpret: the effect of the EM force 
is such that the motion of e and p is tightly coupled to- 
gether, so that their effective mass per particle is m p /2, 
and indeed lighter than m. This results in the diffusion of 
both e and H + against the direction of the gravitational 
acceleration, if helium contribution is negligible. 

We are now ready to include the diffusion of Lithium, 
using already found solutions for V p — Vh and Vh c -Vh. 
The general expression is given by 
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It turns out that, to a good approximation, we can 
neglect the helium contribution, tt,ho/"-h — > 0, thus 
— > 0, m — > m p , and the cumbersome expressions 
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where in the last formula we approximated = 7m 
The direction of the lithium diffusional velocity is far 
from obvious: it depends on the competition of the two 
terms on the r.h.s. of Equation ([8]), and if the friction 
relative to p wins (i. e. small t Lip limit) , the motion of Li + 
ions will trace the motion of ionized fraction of hydrogen 
gas. Indeed, Eq. (fsj) reduces to (Via — Vh)/t p h = — g, 
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or Vl; = V p , if r Li p is the largest parameter. We now 
need additional input with actual size of r^ 1 . 

The scattering of 4 He on p has been calculated in [TT] . 
The value of the transport cross section in the range of 
energies we are interested in, and its weighted average 
over the Maxwellian velocity distribution is given by 
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where T is the temperature of the baryonic fluid, as = 
l/(am e ), is the Bohr radius, and /i 14 = 4m p /5. The cross 
sections of a singly-charged ion on a neutral atom can be 

approximated as a a b — 2.27T ^ """'^" j [T2"] . leading to 
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where a po \(b) is the atomic polarizability of the neutral 
species 6: a po i(H) = %a% and a po i(He) = 1.38a^. Ry 
stands for the hydrogen binding energy, Ry = a 2 m e /2 ~ 
13.6 eV. We should note that p-H scattering is in practice 
a more complicated process due to the identical nature of 
the nuclei involved, and a far more elaborate treatment of 
the p-H cross section can be found in [15] . However, for 
the accuracy of our discussion, we shall still approximate 
it with Eq. (10 1. Finally, and most importantly, the 
p-Li + scattering is given by the Rutherford formula, 
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where ^17 = |t71 p , and In A is the Coulomb logarithm. 
For the conditions of primoridal plasma after the recom- 
bination, its value is large, In A ~ 40, and weakly depen- 
dent on temperature. Because of the long-range nature 



of the EM force, Eq. (11) exhibits strong enhancement 
by (Ry/T) 3 / 2 lnA at small velocities/low temperatures. 

We are now ready to determine the sign of the r.h.s. 
bracket in the simplified formula Q: 
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where the abundance of free protons in the primordial 
plasma is the same as the electron ionization fraction 
X e . It is easy to see that for the cosmological parameters 
between recombination and re-ionization, the expression 
( 12 ) is negative. In Figure |l]), we plot the the positivity 
condition on X e -redshift plane, assuming standard rela- 
tions between T, photon temperature T 7 and redshift z. 
The separatrix stays firmly below cosmological X e (T) at 
all redshifts. One can see that even tiny values of X e 
would lead to a tight coupling of lithium to H + , result- 
ing in outward diffusion of lithium. We can also quantify 
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FIG. 1. Post-recombination ionization fraction X e (T) (black 
line) and the separatrix for the direc tion of VLi — Vh relative 
to g. Above the gray line Eq. ( 12 1 holds, and since it stays 



always below the black curve, lithium diffuses "out", leading 
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which arc different only at 0(%) level. Conclusions of 
Eqs. ( 12 ) and ( 13 ) are due to large size of the Ruther- 



ford cross section that overcomes the rarity of H + . Inclu- 
sion of He into this analysis does change the conclusions 
somewhat: while Li + remains tightly bound to H + , the 
outward diffusion of ionized H + is no longer guaranteed. 
We find that for most of the redshift of interest, z > 30, 
both p and Li diffuse out of overdensities not changing 
the qualitative details of the simplified analysis. 

Do we understand the magnitude and sign of possible 
[ 7 Li/H] variations? While we have shown that, quite 
unexpectedly, the direction of lithium diffusion in the 
early Universe after the recombination is against local 
gravitational force, it is clear that it would be difficult 
to create variations in lithium abundance at O(l) level. 
Let us assume to good accuracy that 1 / tli p is the largest 
coefficient, so that motion of Li + and H + are spatially 
linked. Then, using the continuity equations we can tie 
the variation in lithium abundance that develops at red- 
shift Zf to the local halo density, 
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Here p is the average matter density, while phaio is 
the mass density of the halo, which is presumably con- 
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FIG. 2. Level of lithium depletion within primordial dark 
matter haloes at z = 20, with virial gas masses of 10 2 (solid), 
10 3 (dotted) and 10 4 M (dashed). We assume an NFW halo 
|14| with concentration of c = 5 for dark matter, and ignore 
gas overdensity, as expected in the sub- Jeans mass regime. 



tributed to mostly by the cold dark matter in the sub- 
Jeans regime. The Hubble expansion rate, H(z), makes 
appearance in this formula. Notice that while p(z), H(z), 
and T p H depend on time/redshift in a simple calcula- 
ble way, phaio can vary by orders of magnitude, depend- 
ing on the position with respect to dark matter haloes. 
The last line of Eq. ( 15 ) assumes typical overdensity of 
A = Ph a io/p{z) 



200 for collapsed haloes, and ignores 
overdensity of neutral hydrogen, as expected for idealized 
sub- Jeans mass haloes. We have confirmed this estimate 
by direct numerical integration of system Q. 

Judging by these results, one would conclude that the 
post-recombination diffusion-induced variation of lithium 
abundance is going to be small on average. That does not 
mean, however, that such variations will be small every- 
where in the Universe (see Fig. [2]). The central regions 
of the halos, which have the shortest cooling time and 
thus most likely to form the first stars, can have signif- 
icantly higher densities, and thus create 0(1) depletion 
of lithium abundance in a spatially small patch. Fig. ([2| 
may also help explain why the Lithium depletion is larger 
in lowest metallicity stars [3], as they must have formed 
earlier at the centres of more massive haloes (with denser 
cores and thus shorter cooling times), which had experi- 
enced a higher level of depletion. Therefore, in light of 
finding a physical mechanism that can lead to the deple- 
tion of Li in overdensities, we may re-state the lithium 
problem in the following way: in addition to asking "how 
Li got destroyed?", one could also question "how likely 
is that our observations sample a special part of the Uni- 
verse, where Li was depleted before stars formed?" . 

We shall briefly emphasize the limitations of our ap- 
proach: As the halo mass approaches the Jeans mass, 
the gas overdensity will become significant, which in turn 



Supersonic halo mergers lead to turbulence which fur- 
ther suppresses the diffusion. Furthermore, the ioniza- 
tion fractions may differ from the cosmological values 
at significant gas overdensities. On the other end, for 
smaller haloes the supersonic relative velocities of dark 
matter and baryons suppresses diffusion, except within a 
small fraction of cosmic volume Jj>\. Thus, only care- 
ful hydrodynamic simulations of early star formation, 
which include diffusion effects, will be able to confirm 
the real(istic) magnitude of lithium depletion in popula- 
tion II and III stars. 

Finally, we would like to stress that perhaps our 
most interesting and novel result is that lithium remains 
closely tied to ionized fraction of the gas. Therefore, any 
additional forces f Q that act on neutral and charged com- 
ponents differentially, could play a role in creating vari- 
ations in lithium abundance. Interesting candidates for 
creating such forces are radiation pressure/stellar winds 
from first stars, and possibly primordial magnetic fields. 
This whole scope of issues deserves close attention due 
to the continuing interest in the lithium problem. 

We would like to thank Dr. J. Pradler for helpful dis- 
cussions. 



suppresses t p h m Eq. (151, and thus the depletion rate. 
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